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Abstract 

We examine of the capability of the Next Linear Cohider to determine the mass as 
weh as the couphngs to leptons and 6-quarks of a new neutral gauge boson, Z', below 
direct production threshold. By using simulated data collected at several different 
values of y/s, we demonstrate how this can be done in a model-independent manner 
via an anonymously case approach. The importance of beam polarization to the success 
of this program is discussed. The procedure is shown to be easily extended to the case 
of top and charm quark couplings. 
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1 Introduction 



While the Standard Model(SM) is in full agreement with all experimental data[^, it is 
generally believed that new physics must exist at a scale not far beyond the reach of existing 
accelerators. Associated with this scale may be a host of new and exotic particles. A new 
neutral gauge boson, Z', is the most well-studied of all exotic particles and is the hallmark 
signature for extensions of the SM gauge group. Current direct searches for the existence of 
such particles at the Tevatronp| suggest that their masses must be in excess of 500-700 GeV 
depending upon their couplings to the SM fermions and their kinematically accessible decay 
modes. If such a particle is found at future colliders the next step will be to ascertain its 
couplings to all of the conventional fermions. In this way, we may hope to identify whether 
this new particle corresponds to any one of the many Z"s proposed in the literature or 
is something else entirely. At hadron colliders, a rather long list of observables has been 
proposed over the years to probe these couplings - each with its own limitations^, It 
has been shown under idealized conditions, at least within the context of E'e-inspired models, 
that the LIIC(y^ = 14 TeV, 100/6^^) will be able to extract useful information on all of the 
Z' couplings for Mz' below ~ 1 — 1.5 TeV. It is not clear, however, how much of this program 
can be carried out using realistic detectors at the LHC|^ and how well it generalizes to other 
extended gauge models since detailed simulation studies have yet to be performed. 

At the NLC, when ^/s < Mz' (the most likely scenario for a first generation machine 
given the Tevatron bounds) a Z' can only manifest itself indirectly as deviations in, e.g., cross 
sections and asymmetries from their SM expectations. This is analogous to the observation 
of the SM Z at PEP/PETRA/TRISTAN energies through deviations from the expectations 
of QED. Fortunately the list of useful precision measurements that can be performed at the 
NLC is reasonably long and the expected large beam polarization(P) plays an important 
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role-essentially doubling the number of useful observables. In the past, analyses of the ability 
of the NLC to extract Z' coupling information in this situation have taken for granted that 
the value of M^/ is already known from elsewhere, e.g.^ the LHC[^ |^. (In fact, one might 
argue that if a 1 TeV Z' is discovered at the LHC, a future lepton linear collider designed 
to sit on this Z' must be built and will thus quite easily determine all of the Z' couplings in 
analogy to SLC and LEP.) Here we address the more complex issue of whether it is possible 
for the NLC to obtain information on couplings of the Z' if the mass were for some reason a 
'priori unknown. In this case we would not only want to determine couplings but the Z' mass 
as well. We will limit our discussion below to the e^e~ channel and ignore the additional 
information available through e~e~ collisionsP|. 

If the Z' mass were unknown it would appear that the traditional NLC Z' coupling 
analyses would become problematic. Given a set of data at a fixed value of i/i which shows 
deviations from the SM, one would not be able to simultaneously extract the value of M^/ 
as well as the corresponding couplings. The reason is clear: to leading order in s/M|,, 
rescaling all of the couplings and the value of Mz' by an overall common factor would leave 
the observed deviations from the SM invariant. In this approximation, the Z' exchange 
appears only as a simple contact interaction. Thus as long as ^/s < Mz', the only potential 
solution to this problem lies in obtaining data on deviations from the SM at several, distinct 
a/s values and combining them into a single fit. It is clear from the beginning that all of the 
set of ^/s values chosen for this analysis cannot lie too far below the Z' mass otherwise we 
would always remain in the contact interaction limit. It must be, for at least one of the ^/s 
choices, that sub- leading terms of relative order s/M"^, are of numerical importance. This 
suggests that the maximum value of ^/s within this set should only be about a factor of 2-3 
lower than the Z' mass. 
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Here we report on the first analysis of this kind, focussing on observables involving 
only leptons and/or 6-quarks. In performing such an analysis, we need to know how many 
^/s values are needed. We need to know how we distribute the integrated luminosity(£) to 
optimize the results. Similarly, we must address whether such an analysis can be performed 
while maintaining model-independence. In this initial study we begin to address these and 
some related issues. It is clear that that there is a lot more to be done before we have the 
answers to all these questions. 



2 Analysis 

In order to proceed with this benchmark study, we will make a number of simplifying as- 
sumptions and parameter choices. These can be modified at a later stage to see how they 



influence our results. (The basic analysis follows that discussed in 0.) In this analysis 
we consider the following ten observables: the total production cross sections for leptons and 
6-quarks, cr^^fe, the corresponding forward-backward asymmetries, A^p^, the left-right asym- 
metry obtained from flipping the initial electron beam polarization, A^^, and the polarized 
forward-backward asymmetry, Ap^{i,b). For r+r~ final states, we include the average r 
polarization, < >, as well as the forward-backward asymmetry in the r polarization, 
Pf^. Other inputs and assumptions are summarized as follows: 

e,/x,r universality ISR with y/^/y/s > 0.7 

P = 90%, SP/P = 0.3% SC/C = 0.25% 

eb = 50%, Ub = 100% 1^1 > 10° 

ee,/^,T(cr) = 100%, eri^Pr) = 50% Neglect t-channel exchange in e~^e~ e+e~ 



Of special note on this list of assumptions are: (i) a 6-tagging efficiency(et) of 50% for a 
purity(n6) of 100%, (ii) the efficiency for identifying all leptons is assumed to be 100%, 
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although only 50% of r decays are assumed to be polarization analyzed, {iii) a 10° angle cut 
has been applied to all final state fermions to mimic the anticipated detector acceptance, 
(iv) a strong energy cut to remove events with an excess of initial state radiation(ISR) has 
been made-this is critical since events with lower effective values of ^/s substantially dilute 
our sensitivity, (v) it has been assumed that both the beam polarization(P) and machine 
luminosity(£) are both well measured. It is important to note that we have not included the 
t-channel contributions to e^e~ e~^e~ in these calculations. In addition to the above, final 
state QED as well as QCD corrections are included, the 6-quark and r masses have been 
neglected, and the possibility of any sizeable Z — Z' mixing has also been neglected; this 
is an excellent approximation for the Z' mass range of interest to us given that we are not 
interested in the Z' W^W~ mode. Since our results will generally be statistics limited, 
the role played by the systematic uncertainties associated with the parameter choices above 
will generally be rather minimal, especially in the lepton case. Larger systematics should 
possibly be associated with the 6-quark final states |jll]| but they have been ignored here for 
simplicity. 

To insure model-independence, the values of the Z' couplings, i.e., (f , a)^ 6, as well 
as M^', are chosen randomly and anonymously using a random number generator from 
rather large ranges representative of a number of extended gauge models. Monte Carlo data 
representing the above observables is then generated for several different values of ^/s. At 
this point, the values of the mass and couplings are not 'known' a priori, but will later be 
compared with what is extracted from the Monte Carlo generated event sample. Following 
this approach there is no particular relationship between any of the couplings and there is 
no dependence upon or relationship to any particular Z' model. (We chose to normalize our 
couplings so that for the SM Z, ai = —1/2.) Performing this analysis for a wide range of 
possible mass and coupling choices then shows the power as well as the potential limitations 
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of this technique. 

To get an understanding for how this procedure works in general we will make three 
representative case studies for the Z' mass and couplings, labelled here by I, II and III. 
There is nothing special about these three choices and several other parameter sets have 
been analyzed in comparable detail to show that the results that we display below are rather 
typical. To begin our analysis, let us try choosing three distinct values. Specifically, 
we generate Monte Carlo 'data' at ^/s =0.5, 0.75 and 1 TeV with associated integrated 
luminosities of 70, 100, and 150 /&~^, respectively. These luminosities are only slightly 
larger than the typical one year values as conventionally quoted |jl2| and assumes a reasonable 



time evolution of the collider's center of mass energy. Subsequently, we determine the 5- 
dimensional 95% CL allowed region for the mass and couplings from a simultaneous fit to all 
of the leptonic and 6-quark observables following the input assumptions listed above. This 
5-dimensional region is then projected into a series of 2-dimensional plots which we now 
examine in detail. 



Figs. 1-3 show the results of our analysis for these three case studies compared with 
the expectations of a number of well-known Z' models^]. To be specific we have considered 
(z) the effective rank-5 model(ER5M), which predicts a Z' whose couplings depend on 
a single parameter —-njl < 9 < 7r/2; (m) the Sequential Standard Model(SSM) wherein 
the new W and Z' are just heavy versions of the SM particles (of course, this is not a 
true model in the strict sense but is commonly used as a guide by experimenters); {iii) the 
Un-unified Model(UUM), based on the group SU{2)i x SU(2)g x U{1)y, which has a single 
free parameter 0.24 < < 0.99; {iv) the Left-Right Symmetric Model(LRM), based on the 
group SU(2)l X SU{2)r x U{\)b~Li which also has a free parameter (k = gn/gL > 0.55) 
of order unity which is just the ratio of the gauge couplings and, lastly, {v) the Alternative 
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Figure 1: 95% CL allowed regions for the extracted values of the (a) lepton and (b) 6-quark 
couplings for the Z' of case I compared with the predictions of the model(dotted), the Left-Right 
Model(dashed), and the Un-unified Model(dash-dot), as well as the Sequential SM and Alternative 
LR Models(labeled by 'S' and 'A', respectively.) (c) Extracted Z' mass; only the > branch is 
shown. In all cases the diamond represents the corresponding input values. 
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Figure 2: Same as Fig. 1 but for a different choice of Z' mass and couplings referred to as case II 
in the text. 
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Figure 3: Same as Fig. 1 but for a third choice of Z' mass and couphngs referred to as case III in 
the text. The larger (smaller) allowed region in each case corresponds to P = 80(90)%. 
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Left-Right Model(ALRM), based on the same extended group as the LRM but now arising 
from Eq, wherein the fermion assignments are modified in comparison to the LRM due to 
an ambiguity in how they are embedded in the 27 representation. 

By examining these figures, several things are immediately apparent - the most ob- 
vious being that two distinct allowed regions are obtained from the fit in all three cases. 
This ambiguity is two-fold and not four-fold in that there is a unique choice of 6-couplings 
for a fixed choice of leptonic couplings. (Of course as we might hope, the input values are 
seen to lie nicely inside one of them.) This two-fold ambiguity occurs due to our inability 
to make an absolute determination of the overall sign of one of the Z' couplings, e.g., ag. If 
the sign of were known, only a single allowed region would appear in Figs. 1-3 for both 
leptons and 6-quarks and a unique coupling determination would thus be obtained. Note 
that this same sign ambiguity arises in SLC/LEP data for the SM Z and is only removed 
through the examination of low-energy neutrino scattering. Secondly, we see that the lep- 
tonic couplings are always somewhat better determined than are those of the 6-quark, which 
is due to the fact that the leptonic observables involve only leptonic couplings, while those 
for 6-quarks involve both types. In addition, there is more statistical power available in the 
lepton channels due to the assumption of universality and the fact that the leptonic results 
employ two additional observables related to the final state r polarization. If the 6-quark 
observables had significant systematic errors, these allowed regions would become larger still. 
Thirdly, we see in particular from Figs, la-b and 3a-b the importance in obtaining coupling 
information for a number of different fermion species. If only the Fig. la(3b) results were 
available, one might draw the hasty conclusion that an E'g-type Z' had been found. Fig. 
lb(3a) clearly shows us that this is not the case. Evidently none of the Z"s associated with 
cases I-III correspond to any well-known model. Fourthly, we note that changing the beam 
polarization from 90% to 80% does not appreciably alter our results. Lastly, as promised. 
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the Z' mass is determined in all three cases, although with somewhat smaller uncertainties 
in case II. It is important for the reader to realize that there is nothing special about any 
of these three particular cases. It is clear from this set of random choices for masses and 
couphngs that this procedure should be viable for Z' masses up to about 2 TeV for the set 
of integrated luminosities that we have chosen unless both of the leptonic Z' couplings are 
accidentally small resulting in a reduced sensitivity to the existence of the Z' . 

It is straightforward to extend this analysis to fermion final states other than lep- 
tons and the 6-quark. The extrapolation to charm is the most obvious in that apart from 
identification efficiencies and potentially larger systematic errors there is little difference 
in performing the five-dimensional fit with either ft or c since the fermion couplings were 
randomly chosen. (Of course we might imagine, however, now doing a more ambitious 
seven-dimensional fit with all of the couplings being allowed to float.) The extension to 
t-quarks would be less straightforward due to their large mass and rapid decay to hW . In 
principle, however, the same set of observables could be constructed for top as was used in 
the 6-quark analysis above. 

Of course, the clever reader must now be asking the question 'why did we start off 
using 3 different values of \fs - why not 2 or 5?' This is a very important issue which we can 
only begin to address here. Let us return to the mass and couplings of case I and generate 
Monte Carlo 'data' for only two values of i/s=0.5 and 1 TeV with luminosities of L— 100 
and 220 respectively, thus keeping the total L the same as in the discussion above. 

Repeating our analysis we then arrive at the '2-point' fit as shown in Fig. 4a; unhke Fig. la, 
the allowed region in the leptonic coupling plane does not close and extends outward to ever 
larger values of v^, a^; we find that a similar result occurs for the 6-quark couplings which are 
even more poorly determined. The corresponding Z' mass contour is also found not to close, 
again extending outwards to ever larger M^/ values. We realize immediately that this is just 
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Figure 4: Failure of the method in case I when data is taken at (a) too few ('2-point' fit) or (b) 
too many ('6-point' fit) different center of mass energies for the same total integrated luminosity 
as in Figs. 1-3. The luminosities are distributed as discussed in the text. 

what happens when data at only a single ^/s is available. For our fixed C, distributed as we 
have now done, we see that there is not a sufficient lever arm to simultaneously disentangle 
the Z' mass and couphngs. Of course the reverse situation can also be just as bad. We now 
generate Monte Carlo 'data' for the case I mass and couplings in 100 GeV steps in ^/s over 
the 0.5 to 1 TeV interval with the same total £ as above but now distributed as 30, 30, 
50, 50, 60, and 100 fb~^, respectively. We then arrive at the '6-point' fit shown in Fig. 4b 
which suffers a problem similar to that presented in Fig. 4a. What has happened now is 
that we have spread the fixed C too thinly over too many points for the analysis to work. 
These same results are found to hold for all three cases. This brief study indicates that 
a proper balance is required to simultaneously achieve the desired statistics as well as an 
effective lever arm to obtain the Z' mass and couplings. It is important to remember that 
we have not demonstrated that the '2-point' fit will never work. We note only that it fails 
with our specific fixed luminosity distribution for the masses and couplings associated with 
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cases I-III. It is possible that for 'lucky' combinations of masses and couplings a 2-point fit 
will suffice or it may work if substantially more luminosity is achievable. It is certainly true 
that all cases where at least three values of are used will allow simultaneous mass and 
coupling extraction provided the integrated luminosity is available. Clearly, more work is 
required to further address these issues. 

How do these results change if Mz' were known or if our input assumptions were 
modified? In this case we use as additional input in our analysis the value of Mz' chosen by 
the Monte Carlo and perform four-dimensional fits to the same set of 'data'. Let us return to 
case I and concentrate on the allowed coupling regions corresponding to a choice of negative 
values of ve^; these are expanded to the solid curves shown in Figs. 5a and 5c. (There will 
also be a corresponding region where ve^b are positive, which we ignore for the moment.) The 
large dashed curve in Fig. 5a corresponds to a reduction of the polarization to 80% with 
the same relative error as before. While the allowed region expands the degradation is not 
severe. If the Z' mass were known, the 'large' ellipses shrink to the small ovals in Fig. 5a; 
these are expanded in Fig. 5b. This is clearly a radical reduction in the size of the allowed 
region. We see that when the mass is known, varying the polarization or its uncertainty 
over a reasonable range has very little influence on the resulting size of the allowed regions. 
From Fig. 5c we see that while knowing the Z' mass significantly reduces the size of the 
allowed region for the b couplings, the impact is far less than in the leptonic case for the 
reasons discussed above. Figs. 6a-d show that case I is not special in that similar results 
are seen to hold as well for cases II and III. For both of these cases, as in case I, there is an 
enormous reduction in the size of the allowed region for the leptonic couplings of the Z' but 
the corresponding allowed region for the 6-quark shrinks by about only a factor of two. In 
case III, there is hardly any reduction in the size of the allowed 6-quark coupling region. 
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Figure 5: (a) Expanded lobe(solid) from Fig. la; the dashed curve shows the same result but 
for P = 80%. The smaller ovals, expanded in (b) apply when the Z' mass is known. Here, in 
(b), P = 90(80)% corresponds to the dash-dot (dotted) curve while the case of P = 90% with 
SP/P = 5% corresponds to the square-dotted curve, (c) Expanded lobe(solid) from Fig. lb; the 
dotted curve corresponds to the case when Mz' is known. 
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Figure 6: (a)Expanded lobe(solid) from Fig. 2a; the dashed curve shows the same result but 
for P = 80%. The smaller dotted oval, applies when the Z' mass is known and P = 90%. (b) 
Expanded lobe (solid) from Fig. 2b; the dotted curve corresponds to the case when Mz' is known, 
(c) and (d) show the corresponding results for case III from Figs. 3a and 3b. 
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Figure 7: A comparison of the constraints obtained on the Z' leptonic couphngs in cases 
I-III, in (a)-(c) respectively, both with(sohd) and without (dashed) observables associated 
with beam polarization. In (d) are shown the corresponding 6-quark couplings obtained for 
case II. 
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Just how large a role does large beam polarization play in obtaining our results? This 
becomes a critical issue, especially at higher energies, if our lepton collider is actually a muon 
collider where we may need to trade off luminosity for high beam polarization [1^. As an 
extreme case, we repeat the previous analysis of cases I-III without including the observables 
associated with beam polarization in the fits; luminosities etc remain the same as before. 
The results of this approach are shown in Figs. 7a-d. In case I, shown in Fig. 7a, the loss 
of the polarization dependent observables causes the allowed region in the leptonic coupling 
plane to open up and no closed region is found. Correspondingly, the 6-quark couplings to 
the Z' as well as the Z' mass are constrained but are no longer localized. Somewhat better 
results are obtained in cases II and III, shown in Figs. 7b and 7c. For case II, the size of 
the allowed region essentially doubles in the vi — ai plane and triples in the Vb — ctfe plane, as 
shown in Fig. 7d. The constraints are found to relax in a corresponding manner. In case 
II, while we are hurt by the lack of polarized beams we are still able to carry out the basic 
program of coupling extraction and Z' mass determination - unlike case I. In case III, Fig. 
7c shows that the leptonic couplings are reasonably constrained without beam polarization 
but now neither the Z' mass nor the 6-quark couplings were found to be constrained. From 
these considerations we may conclude that beam polarization is of critical importance to this 
analysis unless the couplings lie in a 'lucky' range. It is clear that for arbitrary values, we 
will not be able to simultaneously obtain mass and coupling determinations without large 
beam polarization. We note, however, that this conclusion can soften dramatically if the Z' 
mass is already known. 

What happens for larger Z' masses or when data at larger values of ^/s becomes 
available? (As stated above, the 'reach' in our coupling determinations was ^ 2 TeV using the 
'data' at 500, 750 and 1000 GeV.) Let us assume that the 'data' from the above three center 
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of mass energies is already existent, with the luminosities as given. We now imagine that the 
NLC increases its center of mass energy to \/s= 1.5 TeV and collects an additional 200 fb"^ 
of integrated luminosity, which corresponds to 1-2 design years. Clearly for Z' masses near 
or below 1.5 TeV our problems are solved since an on-shell Z' can now be produced. Thus 
we shall concern ourselves only with Z' masses in excess of 2 TeV, inaccessible in the lower 
energy study above. Figs. 8a-d show the result of extending our previous procedure-now 
using 4 different ^/s values, for two distinct choices (IV and V) of the Z' mass and couplings. 
These '4-point' results are a combined fit to the data at all four center of mass energies. 
We show the results for both the general case where Mz' is unconstrained as well as when 
it is already determined by other data. (As before, only one of the allowed pair of ellipses 
resulting from the overall sign ambiguity is shown for simplicity.) Note that the Z' input 
masses we have chosen are well in excess of 2 TeV where the LHC may provide only very 
minimal information on the fermion couplings]^. Clearly by using the additional data from 
a run at y^=1.5 TeV this technique can be extended to perform coupling extraction for 
Z' masses in excess of 2.5 TeV. The maximum 'reach' for the type of coupling analysis we 
have performed is not yet determined. It seems likely, based on these initial studies, that 
the extraction of interesting coupling information for Z' masses in excess of 3 TeV may be 
possible for a reasonable range of coupling parameters. 

3 Outlook and Conclusions 

In this paper we have shown that it is possible for the NLC to extract information on the 
Z' couplings to leptons and 6-quarks even when the Z' mass is not a priori known and, in 
fact, determine the Z' mass. This has been demonstrated in a model- independent manner 
by randomly and anonymously choosing the mass and couplings of the Z' and demonstrated 
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Figure 8: Lepton coupling determination for Z"s with masses of (a) 2.33 TeV and (b) 2.51 TeV 
when the mass is unknown(sohd) and known(dottcd) corresponding to cases IV and V discussed 
in the text, (c) and (d) are the corresponding mass determinations which result from the five- 
dimensional fit. These results include an additional 200 fb~^ of luminosity taken at a center of 
mass energy of 1.5 TeV. 
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the power of precision measurements at future linear colliders. The critical step for the 
success of the analysis is to combine the data available from measurements performed at 
several different center of mass energies. For a reasonable distribution of the luminosities 
the specific results we have obtained suggest, but do not prove, that data sets obtained at at 
least 3 different energies are necessary for the procedure to be successful. The mass 'reach' 
for this approach is approximately twice the highest center of mass energy available. Several 
question remain about the optimization of our approach; these will be addressed in future 
work. 
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